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Aqueous acidification is a serious water quality problem facing many regions of 
the world today. Acidification of waters in regions underlain by siliceous bedrock 
may be attributed to the relatively short retention time of waters and relatively 
slow mineral dissolution kinetics, which limit hydrogen ion neutralization within 
the soil environment. Such waters are characteristically low in dissolved salts (dilute) 
and dissolved inorganic carbon (DIC) due to limited contact with the soil environ- 
ment. However, elevated levels of aqueous aluminum have been observed in dilute 
acidified waters [1-3]. Because dilute water systems are inherently low in DIC, 
they are low in inorganic carbon buffering capacity. Consequently other weak 
acid/base systems, such as aluminum and natural organic acids, may be important 
in the pH buffering of dilute waters. 

Important terms used in evaluating acid/base systems are: base neutralizing 
capacity (BNC), which is the equivalent sum of all acids that can be titrated 
with a strong base, and, acid neutralizing capacity (ANC), which is the equivalent 
sum of all bases that can be titrated with a strong acid. Total BNC and ANC 
can be considered as the composite of individual acid/base systems. 

Several researchers have investigated non-inorganic carbon, weak acid/base 
systems in dilute waters. Dickson [4] observed that elevated levels of aluminum 
increased the BNC of Swedish lakes. Waters were strongly buffered by aluminum 
in the pH range 4.5—5.5. The BNC of aluminum was particularly evident when 
acidified lakes were treated with base (limed). Aluminum BNC was comparable 
in magnitude to hydrogen ion and inorganic carbon BNC; therefore, the presence 
of aluminum substantially increased base dose requirements and the cost associated 
with the restoration of acidified lakes. 

Johannessen [5] investigated nonhydrogen/inorganic carbon buffering in Nor- 
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wegian waters. While reiterating the importance of aluminum as a buffer in dilute 
acidified waters, she also evaluated the role of natural organic acids. Natural organic 
matter reduced the degree to which aluminum hydrolyzed in the pH range 5.0— 
5.5, presumably due to complexation reactions and therefore decreased the buffering 
of aluminum. Natural organic matter also participated in proton donor/acceptor 
reactions; the extent to which total organic carbon (TOC) would dissociate/associate 
protons was 7.5 weq-mg organic carbon !. Johannessen [5] concluded that organic 
carbon was the most important weak acid/base system in acidic Norwegian waters 
because of the high organic carbon concentration relative to aluminum. 

Glover and Webb [6] evaluated the acid/base chemistry of surface waters 
in the Tovdal region of southern Norway. The BNC of hydrogen ion was small 
compared to the BNC of weak acid systems. These investigators suggested that 
of the total weak acid BNC, 40-60 ueq-L-! could be attributed to dissolved alumi- 
num and silicon, while 20-50 peq-L~ could be attributed to natural organic acids. 
Solution titrations were characterized as having a major proton dissociation con- 
stant (Ka) of 1078 to 5 X 1077, in addition to some less well defined ionization 
at higher pH values. 

In a comparable study Henriksen and Seip [7] evaluated the strong and 
weak acid content of surface waters in southern Norway and southwestern Scotland. 
In addition to a titrimetric analysis, the aluminum, dissolved silica and TOC content 
of water samples were determined. Weak acid concentrations, determined by a 
Gran [8] calculation, were evaluated by multiple regression analysis. Most of the 
variance in the weak acid concentration could be explained by the aluminum and 
TOC content of the waters. Thus it was concluded that the weak acid content of 
acidified lakes in southern Norway and Scotland was largely a mixture of aluminum 
and natural organic acids. 

Although some research on the weak acid/base content of dilute, acidified 
water systems has been done in Europe, there have been few, if any, studies in 
North America. The intent of this research was to quantify the weak acid/base 
chemistry of dilute acidified lake and stream systems in the Adirondack Mountain 
region of New York state. 


EXPERIMENTAL SITE AND METHODS 


Water samples were collected over an annual cycle at three lakes, Big Moose 
Lake (74°55’ W, 43?30' N), North Lake (74°51' W, 43?43' N) and Little Moose 
Lake (74?55' W, 43?42' N) and eight tributaries, a total of 24 sites. Samples 
were collected approximately every three weeks in the spring, summer and autumn, 
and every four weeks in the winter; a total of 15 sampling dates. 

pH was determined potentiometrically with a glass electrode. Free fluoride 
was measured potentiometrically with a fluoride ion selective electrode. Sulfate 
was measured by the turbidimetric method [9]. TOC was determined by oxidation, 
followed by infrared detection of CO». Basic cations (Ca?*, Mg?*, Na*, K*) were 
determined by atomic absorption spectrophotometry (AAS). Lanthanum was added 
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Table I. Equilibrium Relationships Used in This Study 


Equation Equilibrium Constant Reference 


Hydroxide Ligands 


AJB* + H20 = Al(OH)?* + H+ 1.03 X 1075 13 
AJP?* + 2H50 = AI(OH); + 2H* 7.36 X 1071 13 
AJI?* + 4H50 = AI(OH); + 4H* 6.93 X 10723 13 
Fluoride Ligands 
A]8* + F- = AIF?* 1.05 x 107 14 
AI?* + 2F- = AIFS 5.77 X 1012 14 
A]* + 3F- = AIF; 1.07 X 1017 14 
Alt + 4F- = AIF; 5.37 X 1019 14 
Als+ + SF- = AIF2- 8.33 X 102° 14 
A]9* + 6F- = AIF} 7.49 X 107° 14 
Sulfate Ligands 
APB* + SOZ- = AISO? 1.63 x 10? 14 
Al8* + 2802- = AI(SO4)z 1.29 X 105 14 
Inorganic Carbon Equilibria 
CO; + H20 = H;CO£? 3.39 X 107? 15 
H5CO? = H+ + HCO; 4.47 X 107 15 
HCO; = H* + CO2- 4.68 X 10711 15 


to samples to minimize potential interference by dissolved silica, aluminum, sulfate 
or orthophosphate in the determination of calcium and magnesium [10]. Specific 
conductance was determined using a conductivity bridge [9]. 

Monomeric aluminum was measured using the colorimetric, ferron-ortho- 
phenanthroline method as suggested by Smith [11]. Organic monomeric aluminum 
(Org-Al) was separated using an ion exchange chromatography technique and 
measured as monomeric aluminum [12]. Details of the ion exchange chromatogra- 
phy technique are summarized elsewhere [12]. Inorganic monomeric aluminum 
(Inorg-Al) was calculated as the difference between monomeric aluminum and 
organic monomeric aluminum. 

Water quality data were examined with respect to chemical equilibria. The 
thermodynamic relationships used in this study are summarized in Table I [13— 
15]. The equations used for calculating various aluminum species are summarized 
elsewhere [2,12]. Activity corrections were made using the Debye-Huckel relation- 
ship [15]. 


Titrimetric Determinations 
The ANC of water samples was determined by titration with strong acid. To 


estimate the ANC of a water sample it is necessary to determine analytically the 
equivalence point of the titration. This determination can be accomplished by a 
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Gran plot analysis. The Gran titrimetric technique was developed by Gran [8] 
and has been discussed by many authors [16-19]. The ANC of a monoprotic 
weak acid/base system can be described as follows: 


ANC = [A7] + [OH] — [H+] (1) 


where ANC=sample acid neutralizing capacity (eq-L~4) 
[A^] = weak base content of the water sample (mol-L~) 


As the titration proceeds, the hydrogen ion in the titrant is added to solution 
and the weak base becomes protonated (HA). Eventually the weak base and hydrox- 
ide ion content of the solution will decrease and the hydrogen ion content will 
increase until the ANC is zero ([H*] = [A] + [OH7]). This occurs at the equiva- 
lence point and the volume of titrant required to reach this point in the titration 
is termed the equivalent volume of titrant (Ve). With the equivalent volume of 
titrant it is possible to calculate the ANC of a sample: 


ANC = Ve X Ca X Vo! Q) 
where | Ve— volume of strong acid required to reach the equivalence point of 
M the titration (L) 
Ca = normality of the titrant (eq-L ) 
Vo — original sample volume (L) 


In addition, any point on the titration curve may be defined relative to the 
equivalence point: 


ANC, = (Ve — V) X Ca X (Vo + V) 1 (3) 


where | ANC, — ANC of any point on the titration curve (eq-L^!) 
V — volume of titrant added (L) 


When the volume of titrant (V) exceeds the equivalent volume (Ve), the ANC, 
becomes negative and the following approximation is valid: 


[H+] >> [A] + [OH] (4) 


ANC = —[H*] (5) 


This approximation for ANC (Equation 5) may be substituted into Equation 
3 and rearranged to obtain 


(Vo + V([H*]) = (V — Ve) Ca (6) 
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and finally 
F, = (Vo + V) 1078 = (V — Ve) Ca (7) 
where F, = ANC Gran function 


If F; is plotted as a function of the titrant volume (Figure 1), the linear 
portion of the plot can be extrapolated to Fi = 0, to obtain the equivalent volume 
(Ve). In addition, the slope of this plot should equal the normality of the titrant 
(Ca). 

Two potential sources of error are associated with a Gran analysis. The 
first is termed the weak base error. If in the pH range of linear extrapolation, 
there is a significant concentration of weak base (A^), then the approximation 
made in Equations 4 and 5 is not valid and an error will result in the ANC 
determination. The titration of weak base will result in a low value of F, and a 
decrease in the slope of the linear portion of the Gran plot. In turn, the linear 
extrapolation will produce a lower value of the equivalent volume (Ve) and this 
will result in an underestimate of the measurement of ANC. 


VOLUME OF TITRANT, V 


Figure 1. Schematic illustration of a Gran plot. The linear portion of the curve is extrapo- 
lated to determine the equivalent volume (Ve). The slope of the plot corresponds to the 
normality of the titrant. 
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To illustrate the magnitude of the weak base error, consider a sample volume 
containing a monoprotic weak acid/base system that has a pKa value of 5 and a 
total concentration (CT) of 2 X 1075 mol-L7!: 


CT — [HA] + [A7] (8) 


This sample is titrated with 0.1 N HCl. To evaluate the weak base error of this 
model system, F; was evaluated as a function of titrant volume. The linear extrapola- 
tion of F, was made over three separate pH ranges and apparent values of Ve 
and ANC were calculated (Table II). The analysis was done at pH intervals of 
0.1 pH unit (11 points), and the extrapolation to F, = O was done by a linear 
regression analysis. It is apparent that to minimize the weak base error the pH 
range of linear extrapolation must be considerably below the pKa value of the 
weak acid/base systems present in the sample. It is noteworthy that high correlation 
coefficient (r?) values were observed in the linear regression of the Gran function 
even though there was significant weak base error. It appears that the slope of 
the Gran plot is a better indicator of potential weak base error than the correlation 
coefficient of the linear extrapolation. To minimize weak base error in solutions 
it is best to evaluate the Gran function over a pH range as far below solution 
proton dissociation constants as possible. 

While titrating samples to as low a pH value as possible will minimize the 
weak base error, it will introduce a pH reading error. This is simply due to the 
fact that hydrogen ion activity is greater at lower pH values. A given error in a 
pH reading is magnified at low pH values because the linear extrapolation of the 
Gran plot is made over a greater range and therefore will result in a greater 
error in the ANC determination. 

In this study ANC was measured by incremental strong acid (0.1 N HCI) 
titrations, in equilibrium with atmospheric carbon dioxide, to pH 3.0-3.2. The 
titration curves were analyzed as Gran plots and the ANC was determined. Of 
the weak acid/base systems present in the water samples, natural organic matter 
had proton dissociation constant (pKa) values between 4.0 and 5.0, while aluminum 


Table II. Effect of the pH Range of Linear Extrapolation of the Gran 
Function on Acid Neutralizing Capacity Measurement for a Monoprotic 
Weak Acid/Base System with CT = [HA] + [A7] = 2.0 X 1075 and pKa = 5. 
The Titrant Normality is 0.1 eq-L-!. 


pH Range of Linear Extrapolation 


4.0 — 5.0 3.5 = 4.5 3.0 = 4.0 


p^ 0.9986 0.9999 1.000 
Gran Plot Slope (eq-L~!) 0.0921 0.0987 0.0998 
Measured [A "](eq-L ^?) 11.9 16.2 18.6 


Measured [A ]/Theoretical [A] 0.60 0.81 0.93 
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has a first hydrolysis constant (p*Ks;) of about 5.0. It is evident that there was 
some weak base error and pH reading error in these titrations. However the pH 
range of linear extrapolation, of the Gran plots, was chosen to minimize both 
sources of error. 


ANC Fractionation 


The second set of ANC measurements made in this study involved the fractionation 
of non-inorganic carbon ANC into aluminum and natural organic carbon compo- 
nents. In this analysis three separate sets of titrations were performed on each 
sample. A 100-mL sample was acidified to pH 3 with 0.1 N HCl and purged 
with nitrogen for at least 20 min to remove any inorganic carbon. To the sample 
100 pL of 1 N Tris buffer [tris(hydroxymethyl) aminomethane] was added. The 
Tris buffer served to raise the pH so that strong acid titrations could be made 
and to provide some buffering in the pH 5-7 range so that titration curves could 
be examined in detail. Samples were titrated with 0.1 N HCl using calibrated 
micropipets (20, 50 and 100 pL) to a pH value of approximately 3.0. This procedure 
was done on the original water sample, an aliquot of water sample that was passed 
through an ion exchange chromatography column to remove inorganic monomeric 
aluminum (aluminum desalted) and on a distilled water blank. The ion exchange 
chromatography procedure was the same as that used to separate inorganic and 
organic aluminum forms. 

An example of a fractionated titration curve resulting from this procedure 
is illustrated in Figure 2. The equivalence points of these titration curves were 
evaluated by a Gran analysis [8]. It is evident that there are three distinct titration 
curves. The difference between the total titration curve and the aluminum desalted 
curve was considered to be an estimate of the aluminum ANC. The difference 
between the aluminum desalted curve and the distilled water background curve 
was an estimate of the ANC due to natural organic matter. 

With titration curves for each fraction of aqueous ANC, it was possible to 
use the concept of weak base error to calculate the proton dissociation constant 
(pKa) and the total concentration of proton dissociable/associable material (CT). 
This calculation was made by utilizing a Gran function, corrected for the presence 
of weak base. The Gran function corrected for the presence of a monoprotic weak 
base, with Tris buffer added may be written as follows: 


F; = (VO — V)[H*] = [OH] - [A] — [Tris = Ca (Ve — V) (9) 


where | Ft — weak base corrected Gran function 
Tris” = basic form of Tris (mol-L~*) 


The weak bases (A^, Tris ) may be calculated from equilibrium and mass 
balance relationships: 
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a Che eke 
K+ (H*)/ya- vm 
[Tris-] = — Eres X Krie (11) 


Karis + [H+ YTris 
where y= activity coefficient for an aqueous form of interest 


K — proton dissociation constant of the acid/base form of interest, e.g., 


__ (AWA) 


{HA} 2 


A 


The subscript A^ represents either aluminum or organic matter depending 
on what fraction is being analyzed. Substituting Equations 10 and 11 into Equation 
9, and converting concentrations to activity, the following working Gran relation- 
ship is obtained. 


Fi = (Vo + V)[H*]/yu*] + [OH }/you-] 


(CL. X Ke) E | (CT ris X Kr) ] 
Ka-  (H*]/ya-) (K rris + {H+} Y rris) 


(13) 


The value of Kris is reported to be 8.32 X 10? at 25C [19]; this was verified 
titrimetrically. The titration curve of the aluminum desalted sample was initially 
analyzed in conjunction with Equation 12 to compute K4- and CT,-. This analysis 
required an initial guess of CT4- and Ka-, and the weak base corrected Gran 
function (F$) was calculated. pKa- and then CT4- were incremented until a maxi- 
mum correlation coefficient (r?) was observed. The pKa- and CT,- values at the 
maximum correlation coefficient of the weak base corrected Gran function repre- 
sented the best fit of the titration data to a monoprotic weak acid/base system. 

The values of pK4- and CT,- obtained from the data analysis of the aluminum 
desalted sample titration were designated as pKorg and CTorg respectively and 
used in the modified Gran function: 

"T di fs a3 [_(CTa- X Ka-) 
Fi KVo- V) (H*V ya] + LOH Ver] — [o pos] 


= | (CT tris x K rris) ] _ [ (CTorg x Korg) | (14) 
(Eras X (H*}/yrris) (Korg + {H+} Yoorg) 

The iterative procedure was then applied using Equation 14 to the total 
sample. The values pK4- and CT,- obtained from this analysis were assumed to 
be due to inorganic aluminum. It should be noted that this analysis is not subject 
to the weak base error previously discussed. 


WEAK ACID/BASE SYSTEMS 61 


— Total sample plus tris buffer 


—+- Cation desalted sample plus 
tris buffer 


—0O- Tris buffer background 


Figure 2. A fractionated ti- 
tration curve. The difference 
in ANC between the total 
sample and the aluminum 
desalted sample was consid- 
ered to be the aluminum 
ANC. The difference be- 
tween the aluminum de- 
salted sample and the Tris 
-100 background curve was con- 
sidered to be the organic car- 
bon ANC. 


ACID NEUTRALIZING CAPACITY (weq x1") 


RESULTS 


ANC and pH values in equilibrium with atmospheric CO» obtained from this 
study are represented in Figure 3. Superimposed on this plot is a theoretical titration 
curve of water in equilibrium with atmospheric CO». It is apparent that, in general, 
Adirondack waters had a greater ANC than expected based on equilibrium with 
atmospheric COs. The difference between the measured ANC, and theoretical ANC 
in equilibrium with atmospheric COs, might be considered to be the residual ANC 
(res. ANC), which would represent non-inorganic carbon bases (e.g., aluminum 
or natural organic matter). A plot of the residual ANC as a function of pH is 
presented in Figure 4. The mean residual ANC value was 9.7 ueq-L !. Residual 
ANC was strongly correlated with pH (Figure 3), which might be expected because 
weak acid/base systems would be deprotonated to a greater degree at higher pH 
values. A multiple regression analysis of the TOC, aluminum and dissolved silica 
suggested that TOC, inorganic aluminum and dissolved silica could only slightly 
explain the observed variance (Figure 4). However the results of this analysis 
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Theoretical Inorganic Carbon ANC in Equilibrium 
with Atmospheric CO2 
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Figure 3. ANC values as a function of pH for samples collected in 
this study. Superimposed on this plot is the theoretical inorganic carbon 
ANC in equilibrium with atmospheric carbon dioxide (Pcoz = 1073-5 
atm). The difference between the actual ANC and the theoretical inor- 
ganic carbon ANC at a given pH is considered to be the residual ANC, 
which may be attributed to aluminum or natural organic carbon. 


were not totally unexpected. Only a fraction of the aluminum and organic carbon 
in solution participates in proton consumption reactions, and this fraction is a 
strong function of solution pH. So it is not unreasonable to suggest that multiple 
regression analysis would be a poor means of quantifying residual ANC. 

A more reasonable approach would be to examine the chemical equilibrium 
of the solution with respect to proton donor/acceptor reactions. The monoprotic 
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RESIDUAL ACID NEUTRALIZING CAPACITY (peq xp!) 


-20 
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Figure 4. Residual ANC as a function of pH for samples collected in 
this study. A multiple regression analysis yielded the following empirical 
relationships: Res. ANC = —5.06 X 10-5 + 1.21 X 10-5 (pH), r? = 0.55; 
Res. ANC = —5.52 X 1075 + 1.21 X 10-5 (pH), 7.77 X 107? (TOC); 

r? = 0.57; Res. ANC = —7.02 X 10-5 + 1.44 X 1075 (pH) + 8.68 X 1077 
(TOO) + 1.71 X 10-5 (Inorg-AL, r? = 0.60. 
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weak acid/base systems chosen for this analysis are a simplified representation 
of proton equilibria. Aluminum is a hydrolyzing metal with an insoluble phase 
[13]. Natural organic acids generally exhibit a “smeared” titration curve, which 
has been attributed to several functional groups [21] or electrostatic interactions 
between the same type of functional groups [22]. While the monoprotic weak 
acid/base analysis presented here is a simplification of the actual system, interesting 
information with regard to the nature of acid/base systems in Adirondack waters 
may still be obtained. 


The results of the organic carbon ANC analysis are summarized in Table 
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Table III. Summary of pKorg Values Observed in this Study 


Acidic Streams Acidic Lakes | Nonacidic Lake Total 


Number of Observations 25 12 5 42 
Mean Value 4.44 4.48 4.06 4.41 
95% Confidence Interval +0.41 +0.46 +1.05 +0.55 


* This value is statistically different (P < 0.1) from acidified lake pKorg values. 


III. Proton dissociation constants were comparable with values obtained from other 
investigators [21]. Observed pKorg values for acidified streams were similar in 
magnitude to those of acidified lakes. Limited observations were also made for a 
neutral-pH lake (Little Moose Lake, pH >7); proton dissociation constant values 
were greater (pKa values were lower) than those for acidified lakes (P > 0.1, t 
test). In acidic lakes and streams, aluminum was to a large extent complexed 
with natural organic matter. The increase in observed pKorg values in acidified 
waters may be due to the association of hydrolyzable aluminum with natural organic 
matter. 

Empirical relationships were observed between TOC and CTorg (Table IV). 
Trends for total and stream data sets were statistically significant, while no statisti- 
cally significant trend was observed for lake systems. The poor correlation between 
CTorg and TOC and lake systems may be due to the presence of autochthonous 
organic carbon (e.g., algae) that might not significantly participate in proton dissoci- 
ation reactions. It is interesting to compare the relationship between CTorg and 
TOC observed in this study with values reported by other investigators (Table 
V) [5,7,23]. It is apparent that the value presented by Schnitzer [23] for a “model” 
fulvic acid is considerably greater than values observed for acidified waters. It is 


Table IV. Comparison of CTorg-TOC and CTamInorg-Al Empirical Equations for 
Total, Lake and Stream Data Sets (CT is expressed in mols-L-!, TOC is expressed as 
mg-L~' of C and Inorg-Al is expressed as mg-L~! of AI) 


Standard Error 


I ———— Significance Level, 

Empirical Equation n Slope Intercept ri P 
Total Data 42 

CTorg = 7.54 + 2.65(TOC) 0.34 2.21 0.60 <0.0001 

CT: = 13.7 + 43.0(Inorg-Al) 10.5 2.31 0.30 <0.0002 
Lakes Data 17 

CTorg = 14.1 + 0.42(TOC) 1.04 4.18 0.01 <0.69 

CTa = 8.0 + 52.9(Inorg-Al) 14.0 2.26 0.48 <0.002 
Streams Data 25 

CTorg = 9.40 + 2.57(TOC) 0.42 3.28 0.61 <0.0001 


CTai = 18.8 + 33.3(Inorg-Al) 13.2 3.21 0.22 <0.02 
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Table V. A Comparison of Total, Organic Carbon 
Proton Dissociation Sites Observed in this Study with 
Values Reported by Other Investigators 


Total, Organic Carbon 
Proton Dissociation/Association 


Sites (pu mols-mg Org-C^?) Reference 
20 23 
5.5 7 
RS 5 
2.6 This study 


also noteworthy that the value of total, organic acid proton dissociation/association 
sites per milligram of organic carbon observed in this study is less than values 
observed in studies of acidic lakes in Norway. 

Titration data would suggest (Figure 2) that the aluminum system was dissoci- 
ating at a higher pH value than the organic carbon system. The mean pKa value 
of the inorganic aluminum system (the total less aluminum desalted fraction) was 
5.55 [P(4.82 < p < 6.28) = 0.05]. This agreed well with AlCl; solutions, which 
were also titrated and analyzed as a monoprotic system (pKa: = 5.45). Empirical 
relationships were observed between CT4; and the inorganic monomeric aluminum 
content (Inorg-Al) of water samples in total, lakes and streams data sets (Table 
IV). 

While fitting aluminum titrations to a monoprotic dissociation model may 
yield interesting empirical observations, it is not a realistic model due to the fact 
that aluminum is a hydrolyzing metal. Aluminum chemistry was also evaluated 
using water quality observations from lake and stream systems. Aquo aluminum 
levels were calculated using inorganic monomeric aluminum, pH, sulfate and free 
fluoride values as well as pertinent equilibrium constants (Table I). Aquo aluminum 
levels closely followed an aluminum trihydroxide solubility model. Values of the 
logarithm of the ion activity quotient (Equation 15) for an aluminum trihydroxide 
model are summarized in Table VI. 


Qr = {Al3+} 
{H*}§ 


(15) 


where Q, = ion activity quotient 
[A1?*] = aquo aluminum activity 


Qp values were very similar for lake and stream systems. These Qp values 
are close in value to the equilibrium constant reported by May et al. [13] for a 
natural gibbsite (log* Kso = 8.77). A plot of pAl** as a function of pH should 
yield a linear relationship with a slope of 3 for an aluminum trihydroxide phase. 
A linear regression of pAl3+ as a function of pH for data obtained in this study 
yields slope values slightly less than three (Table VII). A comparison of the pAI?*- 
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Table VI. Summary of Ion Activity Quotient Values (Q,) for an Aluminum 
Trihydroxide Model (3pH-pA1?*) 


Stream Data Lake Data Total Data 


Number of Observations 116 205 321 
Mean Value 8.44 8.52 8.49 
95% Confidence Interval S131 3:0.51 +0.88 


pH relationships for lake and stream data sets shows a considerable difference 
between the two systems. The pAl**—pH slope for lake data (2.64) was considerably 
greater than that observed for streams. It is reasonable that lake data would more 
closely resemble an ideal aluminum trihydroxide system than stream data. Streams 
are more subject to dramatic changes in water quality than are lakes. Under high- 
flow, low-pH conditions (spring snowmelt and rainfall events), the dissolution rate 
of Al(OH)3 may not be able to maintain a constant aluminum solubility (Qp). If 
lake surface water samples during spring snowmelt, which were highly undersatu- 
rated with respect to Al(OH)s, are excluded from this analysis, lakewater data 
very closely resemble Al(OH)3 equilibria (Table VII). 

An empirical relationship was observed between the organic monomeric frac- 
tion of aluminum (Org-Al) and TOC (Table VID. The organic monomeric alumi- 
num-TOC relationships were similar for both lake and stream systems. The slopes 
of these two empirical relationships were almost identical for both data sets. How- 
ever, the ordinate intercept of the lakes data set is less than that of the streams 
data set. This implies that lakewater contains more nonaluminum-complexing or- 


Table VII. Comparison of pAI?*—pH and Org-AI-TOC Empirical Equations for Total, 
Lake and Stream Data Sets (Org-Al is expressed in mg-L^! of Al, TOC is expressed in 
mg-L~ of C) 


Standard Error 
————— Significance Level, 
Empirical Equation n Slope Intercept p P 


Total Data 
pAI* = —6.55 + 2.55(pH) 321 0.042 0.229 0.93 «0.0001 
Org-Al = —0.088 + 0.046(TOC) 322 0.0015 0.0086 0.76 <0.0001 
Lakes Data 
pAP* = —6.63 + 2.64(pH) 205 0.042 0.209 0.95 <0.0001 
pAB* = —7.32 + 2.76(pH)* 191 0.027 0.135 0.98 <0.0001 
Org-Al = —0.084 + 0.046(TOC) 205 0.0020 0.0103 0.70 <0.0001 
Streams Data 
pAIB* = —5.50 + 2.40(pH) 116 0.088 0.442 0.88 <0.0001 
Org-Al = —0.046 + 0.046(TOC) 117 0.0022 0.0154 0.77 <0.0001 


a This represents lakewater data with 14 lake meltwater samples excluded. 
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ganic carbon than does stream water, which could be attributed to autochthonous 
organic carbon. 

DISCUSSION OF RESULTS 
The acid/base chemistry of dilute acidified surface waters may be best summarized 


by the system's buffer intensity (£). Buffer intensity is a measure of the system's 
ability to withstand changes in pH: 


d dC 
p= Sa (16) 
dpH dpH 
where B = buffer intensity 


dCs, dC, = quantity of strong base or strong acid (mol-L~) required to 
produce an incremental change in pH of dpH. 


Stumm and Morgan [15] and van Breemen and Wielemaker [24] provide 
excellent discussions on the thermodynamic nature of buffer intensity. The buffer 
intensity of four acid/base systems pertinent to dilute acidified waters [water (H*, 
OH-^), inorganic carbon, aluminum and organic carbon] are: 


Buyo = 2.303 (H*] + [OH ]) (17) 
Brco; = 2.303 ([HCOs] + 4 [CO? D (18) 
Ba = 2.303 (9[AI3*] + 9[AI-F] + 9[A1-SO,] + 4[AKOH)?*] 
+ [Al(OH)3] + [AKOH):) (19) 
Brcoon = 2.303 ([RCOOH] + [RCOO]) (20) 
where [AI-F], [Al-SO,] = sum of the aluminum fluoride and aluminum sul- 


fate complexes, respectively (mol-L-!) 
[RCOOH], [RCOO-]= protonated and unprotonated forms of natural 
organic carbon, respectively (mol-L~*) 
The development of these equations is summarized by van Breemen and 
Wielemaker [24]. To apply these buffer intensity relationships to Adirondack surface 
waters the following conditions were assumed: 


1. The system was in equilibrium with atmospheric carbon dioxide (Pcoz = 
10735 atm). 

2. The system had a fluoride content of 0.1 mg-L~! and a sulfate content of 
6.0 mg-L71. 

3. The system was in equilibrium with an aluminum trihydroxide phase with 


a solubility constant of *Kso = 10949 = {Al3*}/{H*}?. 
4. The system has an ionic strength of 3 X 1074. 
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5. The system has a TOC content of 5 mg-L~1, which corresponds to CTorg 
= 2.0 X 1075 mol-L^!; this organic acid system has a monoprotic proton 
dissociation constant of Korg = 107** = (H*)(RCOO-]/[RCOOH]. 


The above conditions were considered to be representative of the mean water 
quality data observed in this study and therefore representative of acidic Adirondack 
waters. Appropriate thermodynamic relationships were applied to Equations 17— 
20 and individual buffer intensity values were computed for a range of pH values. 
A negative logarithm of the buffer intensity (p 8) diagram is presented in Figure 
5. At low pH values (pH <5), the system is strongly buffered by the aluminum. 
As the pH increases, aluminum buffer intensity decreases dramatically. In the 
pH range 5-6 the buffer intensity reaches a minimum. The buffer intensity of 
water ([H*]) and organic solutes are significant in this pH range. As the solution 
pH value increases above 6 the inorganic carbon system becomes the predominant 
buffering system. These calculations were based on the assumption that the system 
was in equilibrium with atmospheric carbon dioxide (Pcoz = 10-55 atm). As a 
result, the inorganic carbon buffering illustrated here would represent the maximum 
that could be observed. If the system was not in equilibrium with atmospheric 


Aluminum 


Inorganic 
carbon Figure 5. Buffer intensity dia- 
gram for dilute Adirondack wa- 
ter systems. Equilibrium with 
aluminum trihydroxide 
(pKso — 8.49), organic solutes 
Organic (CTorg = 2 X 1075, 

solutes pKorg — 4.4) and atmospheric 
carbon dioxide (Pcos = 1073.5 
atm) was assumed. The dashed 
line away from pKorg suggests 
that this model is not entirely 
representative of proton equilib- 
ria of natural organic carbon. 
Multiple functional group or 
electrostatic interaction models 
may more realistically simulate 
proton dissociation away from 
pH pKorg. 
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CO; the inorganic carbon buffering intensity would be smaller in magnitude. These 
calculations would suggest that the lower limit of pH values observed in acidic 
Adirondack waters was regulated by aluminum buffering. 


Significance and Application with Respect to 
Base Addition of Acidified Waters 


Acidification of soft waters is a serious water quality problem. The results of 
this investigation may be applied to an understanding and possible solution of 
the problem of surface water acidification. One such solution might be the restora- 
tion of waters by base addition. The addition of base (CaO, Ca(OH)2, CaCO3) 
for water quality restoration has been practiced in regions experiencing acidification 
[24,25,26,27]. In order to achieve adequate neutralization it is necessary to consume 
all the system base neutralizing capacity (BNC), including aluminum and organic 
carbon BNC. To illustrate this consider a “typical” Adirondack lake (Table VIII) 
whose composition is based on mean observations obtained from this study. Utiliz- 
ing this “typical” lake, estimates of the amount of base required to neutralize 
the system to pH 6.5 can be made, and are presented in Table IX. The calculations 
presented in Table IX were made by assuming: 


1. The aqueous system was in equilibrium with atmospheric carbon dioxide 
(Pcos = 10735 atm) 

2. The 5 mg C X L~! corresponds to 2 X 1075 mols X L™ and the organic 
carbon system has a pKa — 4.4. 

3. Interactions with the sediment are negligible. 


It is apparent (Table IX) that the aluminum and organic carbon systems 
represent a significant fraction (~50%) of the aqueous BNC of Adirondack lake 


Table VIII. Typical Adirondack Lakewater Composition 
(mg-L™ of Al, Except pH and TOC) with Respect to 
Parameters that Exert Base Neutralizing Capacity 


Parameter Value 
pH 4.95 

Inorganic Monomeric Aluminum 0.200 
Aluminum Fluoride Forms 0.105 
Aluminum Sulfate Forms 0.005 
Aquoaluminum 0.040 
Al(OH)?* 0.030 
Al(OH); 0.020 
TOC 5.08 


a mg-L'! C. 
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Table IX. Amount of Base Required to Neutralize BNC 
Components of Typical Adirondack Lakewater to pH 6.5 


Base Required Percent of Total 
BNC Component (10-5 eq-L^?) Base Required 
Hydrogen Ion 1.1 23 
Inorganic Carbon L3 27 
Aluminum 2.0 42 
Organic Carbon 0.4 8 
Total 4.8 100 


systems. Therefore aluminum and organic carbon BNC should be considered when 
evaluating base addition requirements of acidified lakes. 


It is noteworthy that Scheider et al. [26] suggest that sediments represent a 


significant base sink in acidic lakes. This study did not address this apparently 
important component of BNC. 


CONCLUSIONS 


Based on this study it may be concluded that: 


É 


Natural organic carbon and inorganic aluminum are significant weak acid/ 
base systems in dilute acidified waters of the Adirondack region of New 
York state. 

Natural organic acids were fit to a monoprotic proton dissociation constant 
model (pKa = 4.41) and the total, organic carbon proton dissociation/associa- 
tion sites were observed to be empirically correlated to total organic carbon 
concentration. 

Aquo aluminum levels appear to fit an aluminum trihydroxide solubility 
model. 

Aluminum buffering appears to control the lower limit of pH value observed 
in Adirondack waters. 
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